The Two-impurity Anderson Model Revisited: Competition between Kondo Effect and 
Reservoir-mediated Superexchange in Double Quantum Dots 
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We study a series-coupled double quantum dot in the Kondo regime modeled by the two-impurity Anderson 
model and find a new conduction-band mediated superexchange interaction that competes with Kondo physics 
in the strong Coulomb interaction limit. Our numerical renormalization group results, complemented with the 
higher-order Rayleigh-Schrodinger perturbation theory, show that the novel exchange mechanism leads to clear 
experimental consequences that can be checked in transport measurements through double quantum dots. 
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Quantum dots (QDs) Jlj] behave as quantum impurities, 
thus they exhibit the Kondo effect, the most spectacular man- 
ifestation of which is the transition from near-zero conduc- 
tance due to Coulomb blockade to perfect transmission [0, H] 
as the temperature is lowered well below the Kondo tempera- 
ture, Tk ■ More complex setups allow tailored realizations of 
strongly correlated electron systems in new situations. Dou- 
ble quantum dots (DQDs) 0. B S 01, for instance, are a 
minimal system for studying a lattice of magnetic impuri- 
ties in a tunable environment. In the two-impurity Kondo 
model, the competition between the Kondo effect and the anti- 
ferromagnetic (AF) interaction between impurities leads to a 
second-order quantum phase transition (QPT) |8]. It is also 
known that the QPT is replaced by a crossover if electrons can 
tunnel between the impurities [9[]. In series-coupled DQDs, 
the most important contribution to the total AF exchange cou- 
pling is the superexchange coupling Jtj w At 2 /U, where t 
is the inter-dot hopping and U the on-site interaction. This 
superexchange vanishes^ Jtj — ► 0) as U — > oo, which leads 
to the common belief |4|, |7[] that the Kondo physics always 
prevail in the large- U limit. In this Letter, we present numer- 
ical renormalization group (NRG) calculations which confute 
the above claim and unambiguously demonstrate the presence 
of a new type of spin-spin interaction which competes with 
Kondo physics in the large- U limit. Our main findings are 
illustrated in Fig. 0] where the peak position of the linear con- 
ductance through a DQD versus t exhibits a striking and unex- 
pected dependence on V, the QD level broadening due to the 
coupling to the reservoir. Interestingly, the novel spin-spin in- 
teraction reported here leads to observable experimental con- 
sequences like zero conductance G w in regions where uni- 
tary transport G ~ 2e 2 /h is expected and vice versa. Intu- 
itively, we can understand this spin-spin interaction in terms of 
virtual tunneling events involving conduction-band electrons; 
see Fig.[T] The strength of this spin-spin interaction, which we 
call conduction-band-mediated superexchange, Jj, can be es- 
timated in the large-?/ limit using Rayleigh-Schrodinger per- 
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FIG. 1: The sequence of six virtual processes of the most signifi- 
cant contribution to the conduction-band-mediated exchange Ji in 
the large-f/ limit. Here the circles (rectangles) denote dots (leads). 



turbation theory (RSPT) as 

4< 2 r 2 rD 
Ji = — 5— 
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for t -C |e| where e < is the single-particle energy level of 
QD. Here Ef=0 and D are the Fermi energy and the half- 
width of the conduction band, respectively. For a wide band 
(|e| < D), it simplifies to Jj « ct 2 T 2 /\e\ 3 with a constant 
c = 8(1 - m2)/37r 2 « 0.083. Interestingly, a single lead is 
sufficient to induce a finite exchange Jj ^ 0, as seen in Fig.Q] 
Remarkably, this high-order tunneling process is able to affect 
the transport properties of DQDs. In what follows, we present 
in detail our NRG study and discuss thoroughly the conditions 
under which Jj can compete with Kondo physics in regions 
where Jtj vanishes. 

Model and methods. — We model the DQD by the two- 
impurity Anderson model: 



(2) 



Here destroys a spin-/i electron with energy et in the 
lead I = 1,2, and dgn an electron in the dot I = 1,2; ni = 
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di^dgfj, is the occupation of the dot I, For simplicity, we 
assume reflection symmetry with respect to the interchange 
1^2 of leads and dots. The single-particle energy e on each 
dot is tunable by gate voltages. The hybridization between the 
dot and lead is characterized by V = irpV 2 , with a flat-band 
density of states p = 1/2D. Throughout this work, we fix 
e = — 0.1D and focus on the Kondo regime with localized 
level, r |e| and large charging energy U > 2\e\. 

We solve the Hamiltonian with the standard NRG proce- 
dure [Til [l2tl . At low temperatures, a Fermi-liquid sys- 
tem is described by an effective Hamiltonian that takes a form 
similar to the original Hamiltonian but with renormalized pa- 
rameters 111311 . Previous works have applied this theoretical 
approach to deal with single dots 111311 . Here we extend it to 
treat the two-impurity problem of a DQD. It is then techni- 
cally convenient to change to the parity basis (even and odd): 
c e {o)kn = (cik/i ± C2k M )/ V2 and d e (o)/i = (dift ± d 2f _ l )/V2. 
In this basis, the Hamiltonian in Eq. (|2) reads as 
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the local spin operator is defined by S s = 
, (T^ft'dl^dgft' with cr being Pauli matrices, and the 
even (odd) level is given by e e( - ) = e =F t. 

The renormalized parameters at the Fermi-liquid fixed point 
(which we denote by appending an asterisk to bare parameter 
symbols) are extracted from the NRG flow diagrams. Ad- 
ditional terms, missing in the bare Hamiltonian, which are 
allowed by the symmetry are generated in the renormaliza- 
tion process. In the absence of the particle-hole (p-h) sym- 
metry, for example, coefficients U s of the term n s ^n s i, hy- 
bridization r s , and energies e s will renormalize differently 
in each parity channel. In the p-h symmetric case we were 
able to determine all interaction terms in the local basis, in 
particular, the exchange coupling J\ 2 between dot spins Si 
and S 2 . Finally, the linear conductance is determined by 
even and odd quasi-particle scattering phase shifts as G — 
(2e 2 /h) sin 2 (<5 e - S a ) 0]. In the following we consider two 
cases: (i) the model with the p-h symmetry and intermediate- 
ly , where the physics is dominated by Jjj and (ii) the large- 
U case where Jjj is suppressed. Analysis of two cases re- 
veals the emergence of a new type of anti-ferromagnetic cou- 
pling J i generated by high-order virtual tunneling events (see 
Fig. [TJ among the electrons localized on dots and the itinerant 
carriers in the reservoirs. 

Particle-hole symmetric, intermediate-U regime. — Let us 
review briefly the main features for the p-h symmetric case 
in Fig. |2] The Friedel-Langreth sum rule imposes the rela- 
tion S e + S = tt for all t and then the average occupation 
per spin channel is one. The physics is governed by the com- 
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FIG. 2: (a) Zero-temperature linear conductance and phase shift dif- 
ference 5 — 5 e — S in the symmetric case (inset) as functions of t/Y, 
U — 2 1 e | . Lines are guidance for eyes, (b) Linear conductance G 
(filled symbols) and phase shift as functions of the scaled exchange 
coupling, (J12 — Jc)T/Ti(t. Results for the phase shift differences 
and for conductance collapse well onto universal curves, the scaling 
functions fitted to the data. 



petition between the Kondo and the AF energy scales. For 
small t, J12 <C Tk and each dot spin forms a Kondo sin- 
glet state with conduction electrons in the neighboring lead, 
while the weak inter-dot coupling yields small conductance. 
Furthermore, the spin-spin correlation (Si • S2) vanishes; see 
Fig.[3ja). As t — > 00, J12 ^> Tk, and the dot spins are locked 
into a local singlet state, thus inhibiting the Kondo effect, and 
consequently the conductance is small in this limit as well. 
As t increases from to 00, the conductance exhibits one pro- 
nounced peak as shown in Fig.[2|a) where the peak position is 
determined by the condition J c (t = t c ) = J i2 ~ Tjc- At this 
point the spin correlation takes the value (Si • S 2 ) = —1/4 
due to the singlet-triplet degeneracypl] . The crossover occurs 
from the Kondo to AF regime J3, Uwhere J12 ~ Tk, and 
the parameters t* and T* become equal at J12 = J c . This 
crossover is known to be well described by a scaling func- 
tion H: S/tt = 0[(Ji2 - J c )T/T K t] with </>(-oo) = and 
4>(oo) = 1 which is in good agreement with our numerical 
results; see Fig.[2jb). 

Large-U regime. — We found that the large-C/ case is gov- 
erned by spin-spin correlations mediated by the conduction- 
band electrons. The linear conductance features a plateau 
starting at t ~ |e| Il4ll and a peak at lower t; see Fig. [4] 
Importantly, the peak reaches the unitary limit and shifts to- 
ward larger t/Y with increasing T as in the symmetric case. 
It must be emphasized that this behavior of the conduction 
peak is wholly unexpected and disagrees with previous ana- 
lyzes for the U — > 00 limit. We therefore continue with a 
more detailed study of this feature. In the large-?/ limit, the 
usual superexchange Jjj oc 1/U vanishes. For this reason, the 
conductance peak has been traditionally ascribed to the for- 
mation of even and odd Kondo states occurring at T* = t*. 
In addition, some theories Ji , L5 ] predicted that T and t are 
renormalized in the same manner, thus the transition from in- 
dividual Kondo states to bonding and anti-bonding states was 



3 



U/D=0.2 




FIG. 3: Zero- temperature properties of the p-h symmetric model (left 
panels) and the p-h asymmetric U = oo model (right panels) for F = 
0.01D. (a,d) Linear conductance and spin-spin correlation (Si -82). 
(b,e) Renormalized parameters, t* , T*, T , |e*|. In the symmetric 
model, r* = r* and e* = 0. (c,f) Renormalized superexchange 
coupling in the local basis (c) and in the parity basis (f). 



expected to happen at T = t. However, the results in Fig.|4]are 
in clear disagreement with this picture: the conductance peak 
is not located at T = t, but rather at t < T and it is shifted 
toward smaller t/T with decreasing T. This suggests the pres- 
ence of some processes which are responsible for the different 
renormalization of T and t. We indeed find an exchange cou- 
pling between the dot spins, which is manifested in the unex- 
pectedly strong anti-ferromagnetic correlation (Si • S2) < 
in Fig. (3|d). At the conductance peaks, this correlation takes 
the value of —1/4. The origin of the exchange coupling at 
U — > 00 is illustrated in Fig. Q] in which the dot spins get to 
be locked into a singlet state by means of virtual tunneling 
processes where conduction-band electrons participate. 

The strength of this conduction-band-mediated superex- 
change coupling Jj in the large-C/ limit is estimated using the 
RSPT; see Eq. (Q]). The RSPT (or projection method) should 
be enough to extract the almost true value of Jj since Jj is due 
to the charge fluctuations, as Jjj is; see Fig.Q] Further scaling 
would affect only the spin fluctuation part related to the Kondo 
effect and not to the local spin interactions. Hence, Ji com- 
puted from the RSPT should be, at least in order of magnitude, 
very close to the one that competes with the Kondo scale. The 
superexchange coupling Jj is generated at the fourth order of 
V in the RSPT scheme and in the limit t -C |e|, U it is of 
the second order in t; see Eq. (HJ. It is remarkable that such 
higher-order processes lead to sufficiently strong spin corre- 
lations to qualitatively affect the transport properties of the 
DQD. The maximal conductance at the crossover region oc- 




FIG. 4: Zero-temperature linear conductance and phase-shift differ- 
ence 5 = S e — So (inset) in the large-t/ case (U = 00). Lines are 
guidance for eyes. The dashed line is obtained from the SBMFT 
calculation @]. 



curs because Jj drives the system from the Kondo phase to 
the anti-ferromagnetic phase at Ji as Tk, in full analogy with 
the symmetric case. This is underpinned by the observation 
that the renormalized e* goes through zero, even and odd hy- 
bridizations become equal, T* = T* = T*, and V* = t*; 
see Fig. [3je). At this particular point, the p-h symmetry is, 
surprisingly, restored at low energies. We also find that the 
magnetic interaction J* a (it is technically difficult to extract 
J*2 unlike in the p-h symmetric case) again exhibits a local 
minimum at the cross over [see Fig. Of)]- The fact that the 
conductance peak is shifted away from T = t can be regarded 
as the experimental proof of the existence of a new type of 
spin-spin correlations not accounted so far. Other theoretical 
techniques like slave-boson decomposition Ji , 15 ] do not de- 
scribe properly these magnetic correlations between the dot 
spins that arise naturally as the interplay between the tunnel- 
ing events and the Coulomb repulsion. In these approaches, 
spin exchange interactions are usually included ad hoc 112I1 . 
Our results suggest that both Jrj and Jj have to be included 
in slave-boson schemes in order to achieve a good agreement 
with NRG calculations. 

Crossover. — Now we address the intermediate regime 
with medium U where both anti-ferromagnetic couplings are 
finite. After a lengthy calculation it is possible to obtain an 
analytical expression for Jj as a function of U flill . The ex- 
pression is too cumbersome to appear here and we just plot 
the final results. The inset in Fig.|5]shows the relative portion 
of J j in the total superexchange versus U. For sufficiently 
large U, Jj surely dominates over Jjj that eventually van- 
ishes as U — > 00. For intermediate values of U < 10D, 
Ji <C Ju the physics is governed by the superexchange in- 
teraction Jjj, while it is Jj that suppresses the Kondo corre- 
lations for larger U. Since Jj is proportional to T 2 , the value 
of U where J/ becomes dominant over Jjj can be lower with 
increasing V. The competition between Jj and Jrj are also 
reflected in the main plot of Fig. [5] where the linear conduc- 
tance peak position versus U is plotted. While for interme- 
diate U the peak position is well identified by the condition 
Ju ~ Tjc (dotted lines), which is in agreement with previ- 
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ferromagnetic interaction for all values of U. Transport exper- 
iments can confirm the presence of high-order superexchange 
coupling by examining the dependence of the peak position 
on the inter-dot hopping amplitude. 
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FIG. 5: The conductance-peak position t c /Y as a function of the 
Coulomb energy U for different values of V: T/D — 0.012 (red) 
and 0.02 (black). The dotted and solid lines are the peak positions 
estimated from the condition Ju — 2.2Tk and Jv + Ji = 2.2Tk, 
respectively. Inset: the relative strength of Ji with respect to the total 
superexchange coupling. 

ous theories, the deviation becomes significant for large U, 
and in the U — ► oo limit the peak position is given by the 
condition Jj ~ Tk- Interestingly, a good fitting for arbi- 
trary value of U is obtained by determining the peak posi- 
tion by imposing the condition J = Jj + Jjj = 2.2Tk, with 
T K = v /rmin(t/,D)/2cxp[27re(l + e/U)/T\. Thus, the 
conduction-band-mediated superexchange explains in a very 
suitable manner (both qualitatively and quantitatively) the sat- 
uration of the conductance peak as U is increased. Our esti- 
mation of the peak position from the RSPT is slightly larger 
than one obtained from the NRG calculations. It is due to our 
neglect of further scaling of Jj with integration of the conduc- 
tion band and the renormalization of Tk due to the inter-dot 
coupling. Both renormalizations should be quite marginal and 
may lead to slightly larger Jj and smaller Tk, with which the 
condition J = 2.2Tr- gives rise to smaller conductance peak 
position. 

Discussion. — The RSPT estimation predicts that Jj be- 
comes comparable to Jjj at U w U c = 4|e| 3 /cr 2 . The ratio 
U c /\e\ is then of the order of a few hundreds in order to retain 
the Kondo effect (|e|/T -C 1). The condition can be satisifed 
through the electric control of the single-particle level in tun- 
able samples, or by making use of ultrasmall samples with 
strong Coulomb interaction. One of the latter canditates is 
a junction made of molecular dimers where U is about a few 
eV. While the interdot coupling of dimers can hardly be tuned, 
the gate voltage control of e and T can nevertheless fulfill the 
crossover condition Jj ~ Tk, allowing to examine the role of 
the superexchange on the transport. 

Conclusion. — We have reported on a new conduction- 
band mediated superexchange interaction, Jj, that competes 
with Kondo physics in the large-[/ limit of the two-impurity 
problem. The novel exchange mechanism brings more unity 
to this problem in terms of magnetic correlations, and the 
transport through the DQD can be analyzed in terms of the 
competition between the Kondo correlation and the anti- 
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